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ABSTRACT 



We have detected emission by the CO 5-4 and 6-5 rotational transitions at 
z = 5.7722 ± 0.0006 from the host galaxy of the SDSS quasar J0927+2001 using 
the Plateau de Bure interferometer. The peak line flux density for the CO 5-4 
line is 0.72 ± 0.09 mJy, with a line FWHM = 610 ± 110 km s'^ The implied 
molecular gas mass is (1.6 ±0.3) X 10^° Mq. We also detect the 90 GHz continuum 
at 0.12 ±0.03 mJy, consistent with a 47K dust spectrum extrapolated from higher 
frequencies. J0927±2001 is the second example of a huge molecular gas reservoir 
within the host galaxy of a quasar within 1 Gyr of the big bang. Observations of 
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J0927+2001 are consistent with a massive starburst coeval with a bright quasar 
phase in the galaxy, suggesting the rapid formation of both a super-massive black 
hole through accretion, and the stellar host spheroid, at a time close to the end 
of cosmic reionization. 

Subject headings: 



1. Introduction 

Understanding the relationship between quasars and their host galaxies has become 
imperative since the discovery of the bulge mass - black hole mass correlation in nearby 
galaxies, a result which implies a fundamental relationship between the formation of super- 
massive black holes (SMBH) and spheroidal galaxies (Ferrarese & Merrit 2000; Magorrian 
et al. 1998; Gebhardt et al. 2000). In our extensive study of the dust and gas content, and 
star formation activity, of the host galaxies of z > 2, luminous (-^-^^^q^ < ~26.1), optically 
selected high redshift quasars (Omont et al. 2003; Carilli et al. 2002; Bertoldi et al. 2003a; 
Petric et al. 2003; Beelen et al. 2006; Wang ct al. 2007a), we have found that roughly 30% 
of these sources are also hyper-luminous infrared galaxies {Lpjn ~ 10^^ L0). This FIR emis- 
sion corresponds to thermal emission from warm dust. Searches for CO emission from these 
sources typically yield molecular gas masses > 10^'' Mq (Solomon & Vanden Bout 2005). 
Such molecular gas represents the fuel for star formation, and can serve as a dynamical tracer 
in a forming galaxy. 

An important finding of our program is that the 30% fraction of FIR luminous quasars 
continues to the highest redshifts, z ^ 6, approaching the end of cosmic reionization, and 
hence pushing toward the very first galaxies and SMBHs (eg. Fan, CariUi, & Keating 2006; 
Fan et al. 2006a). A particularly enlightening example is the highest redshift quasar known, 
J1148-I-5251 at 2; = 6.419. The host galaxy has been detected in thermal dust, non-thermal 
radio continuum, and CO fine emission (Walter et al. 2003; Bertoldi al. 2003b; Carilli et 
al. 2004), with an FIR luminosity (integrated between 42.5;um and 122^m) of 1.3 x 10^^ Lq 
(Beelen et al. 2006), a molecular gas mass of ~ 2 x 10^*^ Mq, and a dust mass of 7 x 10* 
Mq. The dynamical mass to a radius of 2.5 kpc in the host galaxy, as estimated from high 
resolution VLA imaging of the molecular gas distribution, is comparable to the gas mass, but 
much less than the mass predicted by the black hole mass - bulge mass relation measured 
at low redshift (Walter et al. 2004). We have also detected the dominant coohng line of 
interstellar gas, [CII] 158 /xm, from J1148-I-5251 (Maiolino et al. 2005). These observations 
demonstrate that large reservoirs of dust and metal enriched atomic and molecular gas can 
exist in massive galaxies within 1 Gyr of the Big Bang. The current observations suggest 
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active star formation in the host galaxy, with a massive star formation rate of order 10^ Mq 
year~^, although the question remains as to the contribution to dust heating by the AGN 
(Wang et al. 2007a; Li et al. 2007a). 

The source SDSS J092721.82+200123.7 (hereafter J0927+2001) is similar to J1148+5251 
in many ways. J0927+2001 is a quasar selected from the Sloan Digital Sky Survey (York 
et al. 2000), with an optical redshift oi z = 5.77 ± 0.03 (corresponding to tuniv = 1 Gyr, 
assuming a standard concordance cosmology), based on fitting to the (relatively weak) Lya 
+ NV lines, and with Lm = 7 x 10^^ erg s^^ (Fan et al. 2006b; Wang et al. 2007a). The 
black hole mass estimated from both the Eddington limit, and from UV line widths, is ~ 10^ 
M0 (Jiang et al. 2007). 

J0927+2001 has been detected in 250 GHz continuum emission at the IRAM 30m tele- 
scope, with a flux density of 5.0±0.8 mJy (Wang et al. 2006), and at 350 GHz at the Caltech 
Submm Observatory, with a flux density of 18 ±5 mJy (Wang et al. 2007b, in preparation). 
These data, in combination with lower and higher frequency data, reveal a clear FIR excess 
over a standard quasar IR SED. This FIR excess corresponds to thermal emission from warm 
dust with a temperature ~ 47 K. The dust mass is 6.9 x 10^ M0, and the FIR luminosity 
is 1.2 X 10^^ L0. A weak (45 ± 14//Jy) radio continuum counterpart has been detected at 
1.4 GHz, and the rest-frame radio through FIR SED is within the range defined by star 
forming galaxies (Wang et al. 2007a). The observations of J0927+2001 are consistent with 
a massive (stars > 5 Mq) star formation rate of ~ 700 Mq year~^ (Kennicutt 1998), or a 
total star formation rate (stars > 0.1 Mq) a factor 5.6 higher, assuming a Salpeter initial 
mass function, and assuming the warm (47 K) dust is heated by star formation. 

In this letter we present the discovery of CO 5-4 and 6-5 emission from J0927+2001. 
These observations reveal the necessary fuel for star formation, and enable a number of key 
observations in the study of the earliest generation of massive galaxies and SMBH. 

2. Observations and results 

Our observations have been made with the improved Plateau de Bure interferometer 
equipped with the new generation of receivers. With their increased sensitivity (detection of 
sub-mJy lines in a single track), and increased bandwidth (2 of the 4 GHz with the current 
correlator) , these new receivers allow for deep searches for faint redshifted molecular emission 
lines over bandwidths adequate to cover the typical uncertainties of optical redshifts of high 
redshift quasars. 

We searched for CO 5-4 emission from J0927-I-2001 with the Plateau de Bure Interfer- 
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ometer in April and May, 2007. The first track (8 hours) used two polarizations, and two 
frequency settings which covered a total bandwidth of 1.8 GHz (= 5700 km s~^ or A.z — 

0. 13) centered aX z = 5.79 (or 84.87 GHz for redshifted GO 5-4). After a possible line was 
identified in the first observations at 2; ~ 5.77, a second 8 hour track was observed centered 
on the line at 85.12 GHz, covering a total bandwidth of 1 GHz using dual polarization. This 
second track confirmed the 5-4 line with high significance. A third 8 hour track was then 
observed centered on the frequency of the CO 6-5 fine redshifted to 102.10 GHz, which was 
also clearly detected. The spectral resolution for all observations was 2.5 MHz channel"^. 
The phase stability, as measured on the phase calibrator (0851-1-202) was easily adequate for 
coherent integration in the D configuration, for which the interferometric synthesized beam 
had a FWHM ~ 5". Absolute gain cahbration was performed on the star MWC349. 

Figure 1 shows the spectra of the CO 5-4 and 6-5 emission fines from J0927-I-2001 at a 
spectral resolution of 70 km s~^ channel"^. The rms per channel varies sfightly across the 
spectra, but typical values for the 5-4 spectrum are 0.37 mJy beam~^ for 5-4, and 0.52 mJy 
beam~^ for 6-5. 

Both CO emission lines are clearly detected, as well as the underlying continuum emis- 
sion at both frequencies. Gaussian fitting to the lines yield the parameters listed in Table 

1, including: the continuum level, the line peak flux density and FWHM, the velocity in- 
tegrated fine flux, and the line luminosity, L\CO) = 3.25 x 10^'^(1 + z) I Av D\i'^i^^ K km 
s~^ pc^, where lAv is the velocity integrated flux density in Jy km s^'^, is the angular 
diameter distance in Gpc {D^ = 1-23 Gpc aX z = 5.77), and Vobs is the observing frequency 
in GHz (Solomon et al. 1997). A linewidth of 600 ± 70 km s~^, and an LSR line centroid 
of 2; = 5.7722 ± 0.0006, are derived from the spectrum obtained by merging the CO 5-4 and 
6-5 data. 

The measured continuum level for the combined data is 0.12 ± 0.03 mJy. Combined 
with the continuum observations at 250 GHz and 350 GHz, this yields a dust temperature of 
47 K, assuming a modified grey body of index /3 = 1.6 (see analysis in Wang et al. 2007b). 

Figure 2 shows the image of the velocity integrated line emission from the weighted 
sum of the 5-4 and 6-5 fines. The peak of the GO emission is at: (J2000) 09^27^^21. 79" + 
20°01'23.5", or within 0.2" of the optical position. The source appears unresolved at the 
present resolution of FWHM ~ 5" . 
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3. Analysis 

The ratio of the strength of the two CO transitions is: L'^^^j L'^_^ = 1.10 ± 0.27. This 
ratio is consistent with constant brightness temperature (ie. ratio = 1) up to CO 6-5. 
The high CO excitation in J0927+2001 is comparable to that seen in the host galaxy of 
J1148-I-5251 (Bertoldi et al. 2003), and other high redshift quasars (Solomon & Vanden 
Bout 2005; CariUi et al. 2002; Weiss et al. 2007 in prep), as well as in the nuclei of nearby 
nuclear starburst galaxies, such as NGC 253 and Arp 220 (Greve et al. 2006; Bradford et al. 
2003), implying warm (> 50K). dense (> 10^ cm~^) gas. However, given the close spacing of 
the two transitions, observations of a lower order transition (1-0 or 2-1) are needed to better 
constrain the gas excitation conditions. 

We derive the total molecular gas mass (dominated by H2) from the CO luminosity 

using the conversion factor for CO luminosity to total molecular gas mass appropriate for 
Ultra-luminous Infrared galaxies (ULIRGs, or galaxies with Lpjfj > 10^^ Lq year^^; Downes 
& Solomon 1998), or X = 0.8 Mq (K km s^^ pc^)^^, and assuming constant brightness 
temperature from CO 5-4 down to CO 1-0 (ie. constant L'{CO)). The implied gas mass is: 
= 0.8 X L'coi-o = (1.6 ± 0.3) x 10^° Mq. Note that the CO luminosity to gas mass 
conversion factor is a factor five smaller for ULIRGs than it is for normal galaxies, such as 
the Milky Way. For a detailed discussion of the conversion of CO luminosity to gas mass in 
high redshift, FIR luminous galaxies, see Solomon & Vanden Bout (2005). 

The ratio of FIR luminosity to CO luminosity has been used as a metric for the relative 
star formation efficiency in galaxies (ie. proportional to the star formation rate per unit gas 
mass; Gao & Solomon 2004). For J0927+2001, we find: Lpm/L'^o = 650 Lq (K km s"^ 
pc^)"^. In their review of molecular line emission from galaxies at 2; > 1, Solomon & Vanden 
Bout (2005) find that the values of Lfir/L'qq range between roughly 100 and 1000, with a 
mean value of 350 (K km s~^ pc^)~^, for galaxies with typical FIR luminosities between 
10^^ and 10^^ Lq. For comparison, lower luminosity star forming galaxies {Lpj^ < 10^^ Lq) 
have Lpjr/L'qq ratios an order of magnitude smaller (Kennicutt 1998). J0927-I-2001 follows 
this general trend for increasing star formation efficiency with increasing FIR luminosity. 
These results suggest a relatively short gas depletion timescale = (gas mass) / (star formation 
rate) ~ 10^ years, implying a brief, but very intense, starburst. This timescale is comparable 
to the typical lifetime of luminous, high redshift quasars derived from the clustering of 
quasars in the SDSS by Shen et al. (2007). Of course, these arguments assume the FIR 
excess is due to dust heated by star formation, and not the AGN. 

Although the gas mass in J0927-I-2001 is comparable to that in J1148-I-5251, the fine 
width of 600 km is almost a factor two broader than the CO emission from J1148-I-5251, 
and is at the high end of the CO line width distribution observed for high redshift quasar host 
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galaxies (Carilli & Wang 2006; Greve et al. 2005). A simple explanation for a broader line 
would be that the gas disk is more inclined to the sky plane in J0927+2001. High resolution 
imaging of the CO distribution is required to determine the extent, and dynamics, of the 
molecular gas. 

4. Discussion 

J0927+2001 is the second example of an extreme molecular gas mass in a galaxy within 1 
Gyr of the Big Bang. Like Jl 148+5251 at z = 6.42, J0927+2001 shows many characteristics 
of a co-eval massive starburst in the host galaxy of the quasar (Section 1; Wang et al. 
2007b). The observations presented herein reveal the requisite molecular gas reservoir to 
fuel the star formation. As in the case of J 1148+5251, Walter et al. (2004) point out that 
the formation of the heavy elements, and in particular, the ISM processing required to form 
cooler molecular gas and dust, implies that star formation must have commenced very early 
in the host galaxy, at least a few hundred million years prior, or z > 8. 

These results lead to the question: how are such massive galaxies and SMBH formed 
within 1 Gyr of the Big Bang? Li et al. (2007a, b) have addressed this question through 
multi-scale cosmological simulations, including prescriptions for the complex processes of 
star formation and AGN feed-back. They find that early galaxy and SMBH formation is 
possible in rare (comoving density ~ 1 Gpc"^), high density peaks (halo mass ~ 8 x 10^^ 
M0 at 2; ~ 6), in the cosmic density field, through a series of gas-rich, massive mergers 
starting at 2; ~ 14. SMBH formation occurs through a combination of Eddington-limited 
accretion in each progenitor galaxy, plus rapid black hole mergers during galaxy interactions. 
The stellar spheroids are formed in merger-driven, extreme starbursts, with star formation 
rates exceeding 10^ Mq year~^ for short periods (~ 10^ years). The ISM atomic abundances 
rapidly approach Solar in the inner few kpc. Li et al. (2007b) suggest that such systems 
evolve into massive galaxies at the centers of the densest cluster environments seen today 
(~ 10^^ M0). 

Li et al. (2007a) hypothesize that systems such as J1148+5251 and J0927+2001 may 
be in a late starburst phase, where the AGN has recently emerged from its dusty shroud, 
and now potentially dominates the dust heating, even for the warm dust component. The 
question of the relative contribution of star formation and the AGN to the heating of the 
warm dust that produces the observed (rest frame) FIR excesses, will be discussed in detail 
in Wang et al. (2007b), and can be addressed through the observations suggested below. 

One outstanding issue the Li et al. models do not address is the early formation of 
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dust. Such early dust formation remains a puzzle, since the standard ISM dust formation 
mechanism, ie. in the cool winds from evolved low mass (AGB) stars, may require timescales 
longer than the age of the universe at 2; ~ 6. One possible solution is dust formation 
associated with massive star formation (Stratta et al. 2007; Maiolino et al. 2004; Venkatesan, 
Nath, & ShuU 2006; Dwek et al. 2007). 

Such early coeval formation of galaxies and SMBHs has interesting implications on the 
interpretation of the cosmic Stromgren spheres around high redshift quasars, as inferred from 
the Gunn-Peterson absorption spectra (Fan et al. 2006b). Masseli et al. (2007) and Lidz et 
al. (2007) point out that, over the lifetime of the system, the integrated star formation and 
AGN emission may contribute roughly equally to the reionization of the regions immediately 
surrounding the quasar host galaxy. 

The detection of the gas reservoir required to fuel star formation in the J0927+2001 
quasar host galaxy is a key step in studying the formation of massive galaxies at very high 
redshift, but clearly further detailed observations are required to probe the host galaxy and 
SMBH formation in this system. Critical diagnostic observations include: imaging of the 
CO dynamics to constrain the galaxy gravitational mass, and hence the black hole mass - 
bulge mass relation at high z (Shields et al. 2006; Walter et al. 2004), multiwavelength 
imaging of the dust to help constrain dust heating mechanisms, studies of lower order CO 
lines to determine gas excitation and total gas mass (Bertoldi et al. 2004), detection of higher 
density gas tracers, such as HCN, to determine the total mass of gas directly associated with 
star forming clouds (Gao et al. 2007; Krumholz & Thompson 2007), and detection of the 
interstellar gas cooling lines, such as [CII] (Maiolino et al. 2005; lono et al. 2006). The CO 
detection presented herein provides an accurate redshift for the host galaxy as well as direct 
evidence for a highly developed ISM, and hence is the first crucial step for future studies. For 
extreme luminosity objects such as 0927+2001, many of these studies can be performed with 
long integration times on current instruments, such as the Plateau de Bure Interferometer 
and Very Large Array. It will take the Atacama Large Millimeter Array, with its more than 
an order of magnitude increase in sensitivity relative to current millimeter interferometers, 
to perform similar studies on less extreme systems, such as the Ly-a galaxies at 2; ~ 6. 

The Plateau de Bure Interferometer is a facility of IRAM, supported by INSU/CNRS 
(France), MPG (Germany), and IGN (Spain). CC, RW, JW acknowledge support from the 
Max-Planck Society and the Alexander von Humboldt Foundation through the Max-Planck 
Forschungspreise 2005. The National Radio Astronomy Observatory is a facility of the 
National Science Foundation, operated by Associated Universities, Inc. DR acknowledges 
support from the Deutsche Forschungsgemeinschaft (DFG) Priority Program 1177. 
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Fig. 1. — The spectra of the CO 5-4 (upper), 6-5 (middle) and the sum of the 5-4 and 6-5 
(lower) emission from SDSS J0927-I-2001 at 2; = 5.77. Also shown are Gaussian fits to the 
lines with parameters as given in Table 1. All spectra are smoothed to 70 km s~^ channel"^, 
and each independent channel is displayed. The rms per channel varies across the band, 
although typical values in the combined 5-4 and 6-5 spectrum are 0.3 mJy beam~^. 
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Fig. 2. — An image of the sum of the velocity integrated CO 5-4 and CO 6-5 hne emission 
from J0927-F2001 at z = 5.77. Contour levels are: 0.14 Jy km s~^ beam-^ The rms on the 
image is 0.07 Jy km beam~^. The cross indicates the position of the optical QSO. 
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Table 1: Gaussian fitting to the CO line emission from J0927+2001 



Transition 


Continuum 


Line Peak 


FWHM 


lAV 


L'{CO) 




mJy 


mJy 


km s"-*^ 


Jy km s"^ 


K km pc^ 


5-4 


0.12 ±0.03 


0.72 ± 0.09 


610 ± 110 


0.44 ± 0.07 


2.0 ±0.3 X 10^° 


6-5 


0.13 ±0.09 


1.25 ±0.25 


550 ± 150 


0.69 ±0.13 


2.2 ±0.5 X 10^° 



